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ABSTRACT 

An analytical model is developed for the screening of the external magnetic field of a rotating, ax- 
isymmetric neutron star due to the accretion of plasma from a disk. The decrease of the field occurs 
due to the electric current in the infalling plasma. The deposition of this current carrying plasma on 
the star's surface creates an induced magnetic moment with a sign opposite to that of the original 
magnetic dipole. The field decreases independent of whether the star spins-up or spins-down. The 
time-scale for an appreciable decrease (factor of > 100) of the field is found to be ~ 1.6 x 10 7 yr, for 
a mass accretion rate M = lO _9 M0/yr and an initial magnetic moment m = 10 30 Gem 3 which corre- 
sponds to a surface field of 10 12 G if the star's radius is 10 6 cm. The time-scale varies approximately 
as /z 3 8 /M 19 . The decrease of the magnetic field does not have a simple relation to the accreted mass. 
Once the accretion stops the field leaks out on an Ohmic diffusion time scale which is estimated to be 
> 10 9 yr. 

Subject headings: stars: neutron — pulsars: general — — stars: magnetic fields — X-rays: stars 



1. INTRODUCTION 

The decrease of the external magnetic field of accret- 
ing neutron stars has been a long standing puzzle and 
has been explained as being due to Ohmic decay of the 
field (e.g., Goldreich & Reisenegger 1992), crustal mo- 
tion on the star's surface (Ruderman 1991), or "burial" 
or screening of the original magnetic field by the accreted 
matter (Bisnovatyi-Kogan & Romberg 1974). This work 
develops an analytic model for the decrease of the exter- 
nal field due to accretion of plasma. 

For some time after the discovery of pulsars only sin- 
gle stars objects were found. It appeared that pulsars did 
not occur in binaries. Because more that half of all stars 
are in binaries, the occurrence of isolated pulsars was 
explained either by pair disruption during supernova ex- 
plosion leading to pulsar formation, or by the absence of 
SN explosions at the end of evolution of stars in close 
binaries (Trimble & Rees 1971). Bisnovatyi-Kogan and 
Romberg (1974; hereafter BK74) analyzed the evolution 
of X-ray binaries in low-mass systems (e.g., Her X-l) and 
concluded that evolution of such systems should lead to 
the formation of non-accreting neutron star, i.e., radio 
pulsars, in close binary systems. BK74 showed that the 
neutron star rotation is accelerated during the disk ac- 
cretion stage so that a reborn (or recycled) radio pulsar 
should become detectable provided its magnetic field is 
similar to that of isolated pulsars. The absence of radio 
pulsars in close binaries in extensive searches could be 
explained by only one reason (BK74): During the ac- 
cretion stage the magnetic field of the neutron star is 
screened by the inflowing plasma, so that the recycled 
pulsar should have B ~ 10 s - 10 10 G, which is 2 - 4 or- 
ders of magnitude smaller than the field strength of radio 
pulsars. Discovery of the first binary pulsar by Hulse and 
Taylor (1975), and subsequent discovery of more than 50 
recycled pulsars (Bhattacharya & van den Heuvel 1991; 



Lorimer 2001; Lyne et al. 2004) confirmed the conclu- 
sion that recycled pulsars have small magnetic fields as 
predicted by BK74. 

From a statistical analysis of 24 binary radio pulsars 
with nearly circular orbits and low mass companions, 
Van den Heuvel and Bitzaraki (1995) discovered a clear 
correlation between the spin period P and the orbital 
period P or b, as well as between the magnetic field and 
the orbital period. The pulsar period and magnetic field 
strength increase with the orbital period at P or b > 100 
days, and scatters around P ~ 3 ms and B ~ 2 x 10 s 
G for smaller binary periods. These relations strongly 
suggest that an increase in the amount of accreted mass 
leads to a screening of the initial magnetic field, and that 
there is a lowest field strength of about 10 s G. Magnetic 
field screening during accretion has been discussed by a 
number of authors (see, e.g., BK74; Romani 1990; Wijers 
1997; Cheng & Zhang 1998, 2000; Choudhuri & Konar 
2002; Payne & Melatos 2004). The studies by Cheng 
and Zhang and by Payne and Melatos analyze the flow 
and magnetic field evolution within the neutron star in 
contrast with the present work which is concerned mainly 
with the magnetosphere of the star. 

Here we present an analytical model of the screening 
of the magnetic field of a rotating neutron star due to 
the accretion of plasma from a disk. The system is as- 
sumed to be axisymmetric. In the idealized case of a non- 
conducting sphere, the decrease of the external magnetic 
field occurs due to the electric current in the infalling 
plasma. The accretion of this current carrying plasma 
to the star's surface creates an induced magnetic mo- 
ment with a sign opposite to that of the original mag- 
netic dipole. In the more realistic case where the star is 
treated as a conducting sphere, the magnetic field due to 
the accreted current-carrying plasma does not penetrate 
the main volume of the star. As a result the screening 
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Fig. 1. — Three-dimensional view of four field lines going from the star to the Alfven radius va- For the conditions shown the disk 
angular velocity at ta, wk( t a), is larger than the star's angular velocity ui t . 



effect of the plasma is greatly reduced (compared with 
the case of a non-conducting star), and the time-scale 
for appreciable screening is greatly increased. For rep- 
resentative parameters we find that the magnetic field 
of the conducting star decreases significantly on a time- 
scale of the order of 10 7 yr. The screening mechanism 
stops working when the Alfven radius is comparable to 
the radius of the star. 

In the §2 the general picture is described, and the main 
equations are derived. In §3 and §4 the magnetic field 
decrease during plasma accretion is considered for non- 
conducting and conducting spheres, respectively. We cal- 
culate the time-scale for an appreciable decrease of the 
magnetic field of a conducting star which depends on the 
radial thickness of the accreted matter. In §5 we estimate 
the thickness of the accreted matter and the time scale 
for the magnetic field to diffuse out of the star once ac- 
cretion has ceased. In §6 we give the conclusions of this 
work. 

2. THEORY 

We consider disk accretion to a rotating neutron star 
with an aligned dipole magnetic field as indicated in Fig- 
ure 1. That is, we consider an axisymmetric star, mag- 
netic field, and disk. Further, we consider configurations 
which are mirror symmetric about the equatorial plane. 
We use both spherical (R, 9, <p) and cylindrical (r, (f>, z) 
non-rotating coordinate systems. 

For this system the corotation radius is 



density B 2 /8tt (Ghosh & Lamb 1979). This gives 



(GM_\ 



1/3 



1.7 x 10 8 P 2/3 cm , 



(1) 



where, is the angular rotation rate of the star, P = 
2-jt/u* is its period, and the star's mass is considered to 
be 1.4Mq. The Alfven radius r A (or inner radius) of the 
disk is the distance at which the kinetic energy density of 
the matter pv 2 /2 is comparable with the magnetic energy 
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cm 
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Here, fi or ^ 30 = ^/(10 30 G cm 3 ) is the star's mag- 
netic moment, M or Af_ 9 = M/(lCr 9 M Q /yr) is the 
mass accretion rate. The accretion luminosity is L accr — 
GMM/R* w 1.2 x 10 37 M_9 erg/s assuming the star's 
radius i?* = 10 6 cm, while the Eddington luminosity 
for a 1.4M Q star is L Edd w 1.76 x 10 38 erg/s. Thus 
Laccr/ LEdd =~ 0.067Af_g. The dimensionless coeffi- 
cient Ua is of order of one-half and depends weakly on 
the disk parameters (see Ghosh & Lamb 1979; Lovelace, 
Romanova, & Bisnovatyi-Kogan 1995; Long, Romanova, 
& Lovelace 2004). The star's accretion luminosity is 
E a = GMM/R* w 1.2 x 10 37 M_ 9 erg/s with the star's 
radius i?„ assumed to be 10 6 cm. 

The ratio rA/r cr determines the qualitative evolution 
of the system. For < r cr , accretion causes the 
star to spin-up with the rate of increase angular mo- 
mentum dJ/dt w M^GMr^) 1 / 2 , where J — Ioj*, with 
j _ io 45 /45 g cm 2 the star's moment of inertia. On 
the other hand for > r cr , the star spins-down with 
the rate dJ/dt w -M(GMr A ) 1/2 ■ For r cr < r A the star 
spins down but is in the propeller regime where an appre- 
ciable fraction of the accreting matter may be expelled 
from the system by the rotating magnetic field (Illarionov 
& Sunyacv 1975; Lovelace, Romanova, & Bisnovatyi- 
Kogan 1999) or the accretion to the star's surface may 
be highly non-stationary (Romanova et al. 2004) . In the 
propeller regime, the mass accretion rate in the outer disk 
(M in equation 2) may be substantially larger than the 
mass accretion rate to the surface of the star (Lovelace, 
et al. 1999). We consider only cases where > R*, 
with J?* the star's radius. 
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For distances R less than ta, the accreting matter 
moves in a funnel flow which follows approximately the 
dipole magnetic lines as sketched in Figure 2 (see, e.g., 
Romanova et al. 2002). For R <C ta the velocity of 
the matter in the funnel flow is approximately free-fall. 
The cross sectional area of the funnel flow S varies as 
R 3 for a dipole field so that p cx 1/i? 5 / 2 . Consequently, 
pv 2 /2 varies as \jRjl 2 . Thus the kinetic energy density 
of the funnel matter is much less that that of the mag- 
netic field which varies as 1/i? 6 . Therefore, the funnel 
flow is magnetically dominated or force- free for R < ta- 

For the considered axisymmetric system, the magnetic 
field has the form 

B =B p + B 4> 4> , 
with B„ = BrR + BgO . We can write 



B r 



1 



b„ 



1 



i? 2 sin6» 89 ' Rsm.9 dR ' 

where ^(R, 9) = RsinOA^R, 9) is referred to as the 'flux 
function.' Note that ^(R,9) — const is the equation for 
the poloidal projection of a magnetic field line. 

In the force-free limit, the flow speeds are sub-Alfvenic, 
v 2 <C v\ = B 2 /47rp, where va is the Alfvcn velocity. 
In this limit, J x B w 0; therefore, J = AB (Gold & 
Hoyle 1960). Because V-J = 0, (B • V)A = and 
consequently A = A(^>). Thus, Ampere's equation be- 
comes V x B = 47rA(*)B/c. The R and 9 components 
of Ampere's equation imply 



H{-$) = R sin 93^ , and 



dH 4tt , T . 



where H(^f) is another function of Thus, 
H(^) = const along any given field line, and J p = 
{c/4ir)(dH/d^)B p so that (J p • V)iJ = 0. The toroidal 
component of Ampere's equation gives 

A** = -tf(*)^P, (3) 



where 



8 2 1 

+ 



C 2 d 2 



dR 2 ' R 2 d( 2 ' 

and ( = cos 9. This is the Grad-Shafranov equation for 
\E" (see e.g. Lovelace et al. 1986). Note that the lines 
^(R, 9) =const are the poloidal magnetic field lines. Al- 
ternatively the rotation of such a line around the z— axis 
form a flux surface. 
The magnetic moment of the star, 

H= ^- [ d 3 x r x J = /i z 



R<R, 



where 



R 3 dR 



d9 sin 2 i 



</</> 



(4) 



Initially, before there has been appreciable accretion, we 
assume that the magnetic moment ^ = p,(t — 0) is due 
to current flow deep inside the star. Thus the star's ini- 
tial magnetic field is given by ^i(R, 9) = \Xi sin 2 9/R for 
R>R«. Note that A*^ { = 6. 

In the following we show that accretion of current car- 
rying matter to the surface of the star acts to reduce the 
star's magnetic moment. We consider first the case of a 
non-conducting star in §3. In §4 wc consider the realistic 
case of a highly conducting star. 



3. ACCRETION TO A NON-CONDUCTING STAR 

This subsection treats the case of a non-conducting star 
in order to explain our model of the magnctosphcrc. This 
model star has a non-varying dipole field source at its 
center. The magnetic field associated with the accreting 
matter is allowed to penetrate the star's surface R = R*. 

The accreting plasma outside the star is highly con- 
ducting. Consequently, plasma falling onto the star car- 
ries with it an embedded magnetic field and the asso- 
ciated current density as sketched in Figures 2 and 3. 
This plasma infall causes the star's magnetic moment to 
change at the rate 



dfi 7T 3 

dt c 



d9 sin 2 9 [J^vr] 



r»+o 



(5) 



Here, R* + indicates a distance just outside of the star, 
and vr is the plasma infall speed which is approximately 
the free-fall speed vr* = (2GM/i?») 1 / 2 , because the 
star's radius i?» is considered to be significantly smaller 
than the Alfvcn radius ta- This speed is much larger 
than the rotational velocity of the star's surface ~ u>R*, 
because the corotation radius is also considered to be 
significantly larger than R*. We assume that there is a 
strong shock just outside the star's surface which ther- 
malizes the bulk kinetic energy of the flow. 

The total magnetic field \P = fy s + ^ m is the sum of the 
field due to the star, \P S , and that due to current flow in 
the magnetosphcre, For the considered conditions 

where i?» <C ta, we show below that |* TO | <C \& s - We 
have \l/ s = /xsin 2 9/Ria.r from the star, which satisfies the 
A** s = 0. Consequently, equation (3) gives A** m = 
-H(^ s )dH(^/ s )/dty s . Because * m is small, * (without 
the s subscript) in subsequent analysis refers to the field 
of the star. 

The funnel flow as shown in Figure 2 exists close to the 
flux surface ^> = "J/// = h/ta- For a dipole field the fun- 
nel flow hits the star's surface at 9ff w (it^/r^) 1 / 2 . The 
width in \& of the funnel flow is A^fff/^fff Ar^/r^ <C 
1. The radial width is assumed to be of the order of 
the half-thickness of the disk Ha at the Alfven radius 
r a as indicated by MHD simulations of funnel flows by 
Romanova et al. (2002). It is clear that H(^f) is zero 
outside of the funnel flow and that H 2 is a maximum 
in the middle of the funnel flow at — The 
value of is estimated by the condition that B 2 rts B 2 

at the Alfven radius ta- For larger values of (B^/Bp) 2 
at rA, the magnetic field loop would open (Lyndcn-Bcll 
& Boily 1994; Lovelace, Romanova, & Bisnovatyi-Kogan 
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1995). Thus ff 2 t 
show that |* m |/* a « (RJtaY « 1. 
In equation (5) we can use the fact that 



This value can be used to 



J<b = 



c_ _d_H^_ 

4irRsin9 d*T~ 



Thus 



dp 1 2 
dt = l R * 



VR* 



r/2 dH 2 



(6) 



(7) 



R,+0 



where we have combined the contributions of the two 
hemispheres. On the star's surface we have * = 
fism 2 6/R* so that cosfl = (1 - R^/p) 1 / 2 . Thus wc 
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Fig. 2. — Diagram of the poloidal cross-section of the accretion disk and funnel flow. Bp is the poloidal magnetic field, is the toroidal 
magnetic field, J p the poloidal current density, and v p is the poloidal component of the funnel flow velocity. The symbol represents an 
arrow out of the page, while ffi represents the opposite direction. The solid dots indicate a high plasma density. 




V disk 

Fig. 3. — (a) Schematic drawing of an initial dipolc field line (in the x,z plane) and the corresponding twisted field line for the case 
where the angular velocity of the disk at the Alfven radius u>k (ta) is less than the angular velocity of the star u) t . This is opposite to the 
case shown in Figures 1 and 2. The current flow J parallel to the B-field necessary to produce this twist is into the star at high latitudes 
and out of the star at low latitudes. The low-latitude current flow has the dominant screening effect and only this current is shown in (a). 
This current flow along the twisted field line gives rise to a negative toroidal current density near the star as shown in ((b). Panel (b) 
shows the poloidal field component B' arising from the toroidal current J^. The funnel flow is indicated by the small circles and the arrow 
v indicates the flow velocity. 



obtain 



dt 8[i 



dH 2 

(i - R^/fiy/ 2 ~dW 



(8) 



The square root can be Taylor expanded in that 



R*/ta- An integration by parts then gives 



dt 



16^ 2 



/' 

Jo 



(9) 



Thus the star's magnetic moment decreases independent 
of the sign of H and hence independent of the sign of 
^k{ta) - w*. 
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Using the mentioned estimates w /j 2 Jr\ and 

A*// w (h A /r A )^ ff = (h A /r A )(fi/r A ), equation (9) 
gives 

Tt=--^{- A ) {- A ){—) ' (10) 

This equation is not realistic for a neutron star which is 
highly conducting. 

4. ACCRETION TO CONDUCTING STAR 

We now treat the case of a highly conducting star. For 
R < i?* the star is considered to be perfectly conducting 
with the frozen in dipole magnetic field = /Ltj sin 2 9/R. 
The magnetic field associated with the accreting matter 
cannot penetrate the surface R = R*. The shell from 
R = -R* to R a > R* is the layer of accreted matter. 
The layer thickness is R a — i?» = AR a <C R*. The 
accumulated azimuthal current carried by the accretion 
layer is modeled by surface current layer JC<f,(0) at R/c = 
R* + AR a /2. 

The magnetic field in the different radial regions can be 
expanded in terms of solutions of A* "J/ = 0, because the 
contribution of the magnetospheric currents is negligible. 
For the region R > R^ we have 

*i(*,o+ E a «%^' (") 

71=1,3,.. " 



where G n (Q is a Gegenbauer polynomial, C = cos9, and 
the a n are coefficients determined below. The terms with 
even values of n are excluded because they correspond to 
a field which is not mirror symmetric about the z = 
plane. We can take 



G n (0 = (l-O 



2,dPn(0 



(12) 



where P n {C) is the usual Legendre polynomial. The other 
solutions related to the Legendre Q n functions are un- 
physical. One can readily show that the G n obey the 
orthogonality relations 

1 G n {QG m {Q _ 2n(n+l) 



i-C 2 



2n + 1 



(13) 



For example, 



Gx-l-C 2 , 



and 



g 3 = -; 



15 225 



9C 2 



15 



525 



315 



As C 2 -» 1, G n -» n(n + 1)(1 - C 2 )/2. As C -> (the 
equatorial plane), G n — > sin(7m/2)n(n — 2)!!/(n — 1)!!, 
with values of 1, —3/2, and 15/8 for n = 1, 3, and 5. 
For R* < R < Ric, we have 

*- = * i (i?,C)+ E 6 « - j^+t) g «(0 . 

(14) 

so that C) = 0- For R < R*, * = *,(#, £). 

Matching the inner and outer flux functions at i? = i?jc 
gives 

2n+l" 

Rj, 



(15) 



Evaluation of the jump across the R = Rk, surface gives 



V dR 



dR 



E a " 

n=l,3,.. 

^■kRk sin i 



G„(<) 2n + l 



■MO 



(16) 



where /? = R*/R/c < 1, and /C</>(C) is the mentioned 
surface current density. Thus we obtain 

„_4«jgM-^ ^ c Mt)g,(<) (17) 



2n(n + 



(l-C 2 ) 1 / 2 



Note that ai is the contribution to the star's dipole mo- 
ment due to the surface current /Q,, and 03 is related 
to the quadrupole moment due to K.#. Note also that 
1 - (3 2n+1 w (2n + l){AR a /R*)/2 which is justified sub- 
sequently. 

From the considerations of §3, we have 

-s- = (^)«.. 

where J^, is given by equation (6). Thus we find 



(18) 



da n = K +1 vr* 
dt 4n(n + l) 

In particular, 



(1-/3 



2n+lj y 



dP„ dtf 2 



(iai 1 2 , o. 



T/2 (Iff 2 
d# sinfl a 



(19) 



(20) 



which differs from equation (7) by the factor 1 — s . Note 
that we have neglected the gradual variation of the thick- 
ness AR a of the accretion layer. This is justified in §5. 
It is useful in subsequent work to introduce 

A 1 = /j, i + a 1 , An^an/R?- 1 , n = 3,5,.. (21) 

because the A' n s all have the same dimensions as the 
dipole moment. 

4.1. Magnetic Energy 
The magnetic energy of the system is 



W=^ Jd 3 x (B 4 + B') 2 = Wi + W, 

W' = ±- I d 3 x (B') 2 . 
Sir J 



(22) 



Here, B 4 is the initial dipole field and Wi the associated 
magnetic energy which is independent of time. Further, 
B'(i) is the field due to the currents on the surface of 
the star which result from the accretion and W is the 
associated magnetic energy. The cross term involving 
B B' in the first line of the equation vanishes because 
in the volume in which B' is non-zero (R > R*), B^ can 
be written as the gradient of a potential. We can also 
write 

Using equation (16) and equation (13) then gives 



W 



W' = 



En(n + 1) , , . .0 



1 R* 



Rf Ai? a ^ (2n+l)2 
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where 5 n o is the Kronecker delta. We have made the 
approximation 1 - (3 2n+1 w (2n + \)(AR a /R*)/2, which 
is justified below where we estimate AR a . 

Most of the magnetic energy W is in the accretion 
layer on the star's surface. The reason for this is indi- 
cated by the schematic drawing in Figure 3 of the fields 
Bi and B'. The magnetic energy of the B' field outside 
of the star (R> R a ) is 

W'(ext) = \ f 1 d( j~ R 2 dR [(B' R f + [B' e f] 
1 ^ n 2 (n+l) 



Owing to this inequality, 



* n=l,3,.. 



2(2n+ 1) 



(A n - f i i 5 n0 y, (24) 



where we have replaced R a by i?* because AR a <C -R*. 
Retaining only the n = 1 terms gives W'(ext)/W = 
3AR a /(2R*) < 1. 

The 'buried' magnetic field is predominantly parallel to 
the star's surface as indicated by Figure 3. The magnetic 
energy W' of the buried field, denoted B',,, is approxi- 
mately equal to 47ri?^Ai? a (_B||) 2 /87r. Retaining only the 
n = 1 term in equation (23) gives the estimate 

2 R* \A!-im\ 



B' H 



3 AR a i?3 



(25) 



The buried field strength is enhanced over the initial 
dipole field [fa/R 3 ) by a factor w R*/AR a for |Ai| <C /Xj. 



4.2. FieZd at Large Distances from Star 
For large _R, 



i? 



R 3 



In the equatorial plane (( = 0), 

*(iJ,7T/2) = ^ 



2i? 3 



(26) 



(27) 



where ^4_i = fa + ai is the dipole moment at time t. We 
assume here and subsequently that the decrease in the 
dipole moment is limited in the respect that 

3|A 3 |^ 



Ai > 



2r\ 



(28) 



and A 1 > |A„G„(0)|(^/r A )"- 1 for n = 5,7, .. 

Of course, the Alfven radius t"a depends on the mag- 
netic field in the equatorial plane. This field is 

Be (R,n/2) = -^ + 9 -^ + ... (29) 

Owing to inequality (28), \B g (r A , C = 0)| = |-Ai|/r\ to a 
good approximation. For simplicity, we consider that Ta 
depends only on /x or and that the other quantities 
in equation (2) such as M are constant. We then have 



tM 



4/7 



where r A i is the Alfven radius at time t 
(28) can then be rewritten as 

A A. 7/15 / u \ 14/15 

— > 

fa 



A 3 


7/15 


fR*\ 


fa 







(30) 

0. Inequality 
(31) 



r A i \fa 



3/7 



(32) 



We assume that Ka/ta = const <C 1, so that 
A^ff/^ff — h Ail rAi — const <C 1. From §3 we have 
H m = A\lr\. Thus, 

4.3. Dimensionless Variables 

It is natural to measure the magnetic moment A\ in 
units of its initial value fa, 

Ai = — , and A n = — , n = 3, 5, .. (34) 
fa fa 

with the tildes indicating the dimensionless variables. 

We measure * in units of fa/R*, and Br(R*,() m 

units of fa/R 3 . In terms of the dimensionless variables, 

Br(R*,() = -dft/dC,. Further we have 



*ff = ^Ar, and A* ff = ^Ar, (35) 

' Ai 'Ai 'Ai 

from equation (35). 

In subsequent work we drop the tilde's from the A' n s 
and the *'s. 

4.4. Specific Model for H(^>) 

The qualitative dependence of H (vp) is discussed in §3. 
A specific function with this dependence is 



(36) 



We can now rewrite equation (19) in dimensionless form 

as 

cL4„ 2n + 1 J„ ( _ 37 ^ 



dr n(n + 1) ^/7 



Here 
X, 



where = ^(R*,Q, with 

*(o= E ^«g„(o. 

n=l,3,.. 

The dimensionless time variable r is given by 

dt . fr Ai y R* A*// 
7T 



= 4 



(38) 



dr ~\R*J v R * * ff AR a (t) ' 

The initial conditions for equations (37) are Ai(0) = 1 
and A n (0) = for n — 3, 5, ... 

For simplicity, we consider A^fff/^ff to be a fixed 
number small compared with unity. Consequently, the 
solution of equation (37) depends on the dimensionless 
time t and on two dimensionless parameters, £ and e, 

A n = A n (r, £, e) , £ = , e = ^ = , 

rAi * // r A i 



where £<1 and e«l. 



(39) 
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Fig. 4. — Time dependence of the normalized magnetic moment Ai/fii for the special case of only the dipole mode. Here, fj,i is the initial 
magnetic moment of the star, and to is given by equation (43). 



4.5. Single Mode, A\£\ 

For the special case of only one mode A\ ^ 0, equation 
(37) can be approximated analytically. We find 



4/7 



£ 

18/7, 



d( (C 2 -Z)exp[-fc(C 2 -Z) 2 



Here, k = A\' '' / '(2e 2 : £ 2 ) and Z = 1 - £/Ai /7 , where we 
assume that Z > or Ai > £ 7 / 4 . Because e<l, fe 3> 1 
so that equation (40) can be evaluated easily as 



(40) 



2 £ 2 



2 A * /7 ZW ' 



(41) 



Thus wc find 
dAi 

~~dT ~ 



dt 



or A]_ = 1 - - 
16V2 /7 i?» 



> f < 1 
i?* 



(42) 



We have used the approximation Z = 1 which assumes 
Ax » £ 7 /4 « 1. 



Equation (42) can be solved to give 
t = Ml-(l-r) 2( V 7 ] . 

Here, 

28\/2 .R* l 



to 



(43) 



150? v R * e £ 5 <AR a > ' 

is the time-scale for the field decrease. This equation can 
also be written as 

5 



t w 1.5 



-R* r Ai I r A 



vr* h A i \R* J <AR a > 

Here, <AR a > is the time-averaged thickness of the accre- 
tion layer where it is assumed that AR a varies gradually 
with time. For a neutron star of mass M = 1.4M Q and 
radius R* = 10 6 cm, we find 



0.001\ /0.005\ 5 



100 



R*/ <AR a > 



t w 0.78 x 10 5 \ ] | — jr- i | — r7r — | yr 

(44) 

The normalization of £, or equivalently r A i, follows from 
equation (2). The normalization of e, or equivalently 
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hAi/ r Aii is based on the Shakura-Sunyaev disk model 
(1973) which indicates that e depends quite weakly on 
the different parameters (e.g., e oc M 1 / 5 ). The normal- 
ization of R*/ <AR a > is discussed in §5. The shortness 
of this time-scale is due to the fact that only one mode 
is accounted for. 

4.6. Numerical Integrations 

The results of §4.5 for a single mode point up the im- 
portance of using a new time variable in place of t in 
equation (37). That is, we rewrite this equation as 



2n+l 8/7 , 



air 1 n(n + l) 
for n = 1, 3, .., where 



AT{t>) l n {r') 



dr 
dr 



-=A?'\t>). 



(45) 



(46) 



Combining this equation with equation (38) gives 



dt 
d? 



£ 3 v R , <AR a > 1 



R* 



1 13/7 



(47) 



We let ro and Tq denote the values for which A\ < 0.01. 
We have 



to 



f TO dr'Af /7 {r'). (48) 
Jo 



The corresponding actual time when A\ < 0.01 is 



to = 4 



u fe <AR a > 



T (£,e) . 



(49) 



For the case of a single mode (Ai), r' is proportional to 
f (neglecting the weak time dependence of AR a ). For 
this case A\ decreases linearly with f or r'. 

We find equations (45) to be stiff. For this reason a 
predictor-corrector code (Gear 1971) was developed and 
used to solve the equations for 1, 3, and 5 modes. The 
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integrals X n were evaluated using a stretched £ grid with 
10 4 — 10 5 points. Figure 4 shows the numerical solution 
for the case of one mode (Ai). This agrees closely with 
the analytic solution of §4.5. 

Figure 5 shows sample results for cases of three and five 
modes. The time-scales to(£, e) for three and five modes 
are not appreciably different. Combining the results of 
many runs for £ = O.Of — 0.1 and e = 0.02 — 0.1 we find 



to 



f .2 x 10- 4 

C3.6 ,2.7 



(50) 



The integrations become increasingly long as e and £ 
decrease. In view of equation (49) we have 



0.001 



0.005 



100 



"""""t ) {—) U7^?J yr 

(51) 

This time-scale is considerably longer than that for the 
case of the single mode (equation 44). The time-scale 
varies approximately as /z 3 ' 8 /M 19 . Thus for a neutron 
star with an initial surface magnetic field of 2 x 10 12 G 
or jUj = 2 x 10 30 Gem 3 , the time-scale of equation (51) 
increases to 2.2 x I0 8 yr. 

Figure 6 shows the magnetic field lines at a time during 
the evolution and the initial, pure dipole field. Figure 7 
shows the corresponding surface current density on the 
star's surface. 

5. ACCRETION LAYER THICKNESS 

We now give a rough estimate the accretion layer thick- 
ness AR a . We let AM denote the mass accreted after 
time t, which is in a layer of thickness AR a <C -R*. The 
mean density of the layer is thus p ~ AM/(47ri? 2 Ai? a ). 
The change in pressure across the layer in hydrostatic 
equilibrium is p ~ [AM/(47r i? 2 )].g, where g = GMjR\ w 
1.87 x 10 14 cm/s 2 . The equation of state of the neutron 
star matter gives p w k(p) n with k w 2.1 x 10 11 cgs, and 
n w 1.5 (Baym, Bethe, & Pethick 1971). Thus we obtain 



( AM \ ( - n ~ 1 ^ n 



•1.3 x 10 4 



AM 



10- 2 M, 







1/3 



cm 



(52) 



For AR a = 10 4 cm, equation (25) gives an estimate of 
the buried magnetic field of <~ 6.7 x 10 13 G. Thus the 
magnetic pressure in the accretion layer is negligible com- 
pared the matter pressure. 

5.1. Ohmic Diffusion of the Field 

The influence of Ohmic diffusion on the 'burial' of 
magnetic field during accretion to a neutron star has 
been studied by Cumming, Zweibel, and Bildstcn (2001). 
They treat the competition between the inward advec- 
tion of the magnetic field inside the star and the tendency 
of the field to diffusion outward due to the finite conduc- 
tivity of the plasma. For sufficiently large accretion rates 



(such as that considered here) the inward advection is 
found to be larger than the diffusion. 

Magnetic field buried during intensive accretion will 
Ohmically diffuse out of the star after the accretion stops 
(BK74). The time scale for the buried field to diffuse out 
of the star can be estimated as 



tohr, 



(Ai? a ) 2 



V : 



(53) 



where rj is the magnetic diffusivity, and a is the con- 
ductivity. Using Cumming et al.'s §2.2 formula for the 
conductivity and our normalization values we find 



tohm ~ 2.7 x 10' 



/ AR a 



Vio 4 c m y Vio- 2 M e 



AM 



(54) 



Thus it takes a rather long time for the field to diffuse 
out of the star. 

It is possible that the buried magnetic field is subject 
to the interchange or buoyancy instabilities considered 
by Cumming et al. (2001), and by Kato, Fukue, and 
Mincshige (1998), but the detailed analysis is beyond the 
scope of the present work. 

6. CONCLUSIONS 

An analytical model is developed for the screening of 
the external magnetic field of a rotating, magnetized, ax- 
isymmetric neutron star due to the accretion of plasma 
from a disk. The decrease of the field occurs due to 
the deposition of current carrying plasma onto the star's 
surface. This plasma creates an induced magnetic mo- 
ment with a sign opposite to that of the original magnetic 
dipole. The physical mechanism is explained in Figure 3. 
The equations (37) for the field evolution are inherently 
nonlinear and this leads to the generation of higher order 
modes when starting from the lowest order mode. 

The main conclusions from this work are: (1) The field 
decreases independent of whether the star spins-up or 
spins-down (that is, for r cr > or r cr < ta)', (2) The 
time-scale for an appreciable decrease (by a factor > 100) 
of the field is t ~ 1.6 x 10 7 yr for M = 10" 9 M o /yr and 
an initial stellar magnetic moment m = 10 30 Gem and 
it scales approximately as pf 8 /M 19 ; (3) The decrease 
of the magnetic field does not have a simple relation to 
the accreted mass; (4) At late times the magnitude of 
the buried magnetic field is much larger than the initial 
field on the star's surface; and (5) Once the accretion 
stops the field leaks out on an Ohmic diffusion time scale 
which is estimated to be > 10 9 yr. 

The present model has evident limitations which re- 
quire further study: Accreting neutron stars arc not 
expected to have their magnetic and rotational axes 
aligned. The case of small misalignment angles may be 
amenable to analytic treatment. The thickness of the 
accretion layer will in general be a function of position. 

We thank Dr. Dong Lai for a number of valuable dis- 
cussions. We also thank the referee for valuable recom- 
mendations. This work was supported in part by NASA 
grant NAG 5-13220 and NSF grant AST-0307817. 
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